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Abstract 

We determine the fermionic collective modes of a quark-gluon plasma which is anisotropic in 
momentum space. We calculate the fermion self-energy in both the imaginary- and real-time 
formalisms and find that numerically and analytically (for two special cases) there are no unsta- 
ble fermionic modes. In addition we demonstrate that in the hard-loop limit the Kubo-Martin- 
Schwinger condition, which relates the off-diagonal components of the real-time fermion self-energy, 
holds even for the anisotropic, and therefore non-equilibrium, quark-gluon plasma considered here. 
The results obtained here set the stage for the calculation of the non-equilibrium photon production 
rate from an anisotropic quark-gluon plasma. 
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I. INTRODUCTION 



The ultrarelativistic heavy ion cohision experiments ongoing at the Brookhaven Relativistic 
Heavy Ion CoUider (RHIC) and planned at the CERN Large Hadron Colhder (LHC) will study the 
behavior of nuclear matter under extreme conditions. Specifically, these experiments will explore 
the QCD phase diagram at large temperatures and small quark chemical potentials. Based on 
the data currently available from the RHIC collisions there is some evidence that an isotropic and 
thermalized state has been created at times on the order of 1 fm/c [1-4]. The fact that thermal- 
ization proceeds rather rapidly is in contradiction with estimates from leading order equilibrium 
perturbation theory. However, to truly understand how the plasma evolves and thermalizes one 
has to go beyond the equilibrium description and study the dynamics of a non-equilibrium quark- 
gluon plasma. In addition, it is important to know how any deviations from equilibrium affect 
observables so that one might be able to gauge how close the system truly is to being isotropic and 
thermal. 

For example, one would like to know how a momentum-space anisotropy in the distribution 
function of the hard modes would affect observables which are sensitive to the earliest times of 
quark-gluon plasma evolution when the anisotropy is expected to be largest. The best signatures in 
this regard are electromagnetic probes such as photon and dilepton production since these particles 
escape the plasma without strong final state interactions. In order to calculate in-medium photon 
production, however, it is necessary to include the effects of medium-induced fermion masses which 
serve to screen infrared divergences in the calcTilation of production cross sections. In equilibrium 
this can be done self-consistently within the hard thermal loop framework [5-7] and there are now 
many papers dedicated to the calculation of equilibrium photon production at leading and next-to- 
leading order [8-27]. In addition, there have been calculations of electromagnetic signatures from a 
plasma which is not chemically equilibrated [28-36]. However, the problem of photon and dilepton 
production from a quark-gluon plasma which is not isotropic in momentum space has not yet been 
considered. Here we set the stage for such a calculation by computing the quark self-energy in such 
an anisotropic plasma. 

Momentum-space anisotropic distribution functions are relevant because of the rapid longitu- 
dinal expansion of the partonic matter created in a heavy ion collision. This rapid longitudinal 
expansion implies that at proper times r > {pr}^^, where pr is the typical transverse partonic 
momentum of the nuclear wavefunction, the parton distribution functions are oblate in momentum 
space with (px) > {pl)- For RHIC energies this implies that the distribution is oblate for r > 0.2 
fm/c and for LHC r > 0.1 fm/c. Such an anisotropic quark-gluon plasma is qualitatively different 
from an isotropic one since the gluonic collective modes can then be unstable [37-55] . The presence 
of these gluonic instabilities can dramatically influence the system's evolution leading, in partic- 
ular, to its faster isotropization and equilibration. Treating this problem in all of its generality 
is a daunting task. In order to make analytic progress we consider here the limit of high particle 
momentum scale (large pt) and small coupling in order to calculate the fermionic self-energy in 
the hard-loop approximation. 

In two previous papers by Paul Romatschke and one of us [43, 45], we calculated the hard- 
loop gluon polarization tensor in the case that the momentum space anisotropy is obtained from 
an isotropic distribution by the rescaling of one direction in momentum space (corresponding to 
stretching or squeezing of the particle distribution function along a special direction in momentum 
space). In this paper we extend this exploration of the collective modes of an anisotropic quark- 
gluon plasma by studying the quark collective modes using the same framework. Specifically, we 
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derive integral expressions for the quaxk self-energy for arbitrary anisotropy and evaluate these 
numerically using the momentum-space rescaling used in the previous papers. We show for quarks 
there are still only two stable quasiparticle modes and no unstable modes using the momentum- 
space rescaled distribution functions. The result is similar to the case of the fermionic collective 
modes in a two-stream system [56] where it was also found that there were no unstable modes. ^ 

The absence of unstable fermionic modes is expected on physical grounds due to the fact that 
fermion exclusion precludes the condensation of modes; however, it could be possible that, through 
pairing, fermions could circumvent this as has been predicted [57-61] and demonstrated [62] in 
superfluid condensation of cold fermionic atoms. However, this would require a description in 
terms of fermionic bound or composite states which are not included at the level of hard loops so 
we do not expect to find any fermionic condensate-like instabilities using this approximation. This 
is verified via an explicit contour integration of the inverse hard-loop quark propagator for the two 
special cases in which we can obtain analytic expressions for the self-energy. The special cases 
considered analytically are (a) the case when the wave vector of the collective mode is parallel to 
the anisotropy direction with arbitrary oblate anisotropy and (b) for all angles of propagation in 
the limit of an infinitely oblate anisotropy. 

Finally, we present a calculation of the ofi"-diagonal components of the anisotropic fermion self- 
energy using the real-time formalism of quantum field theory. Using this explicit calculation we 
demonstrate that within the hard-loop framework the high-temperature limit of the Kubo-Martin- 
Schwinger (KMS) formula, namely S12 = — S21, holds even for the non-equilibrium configuration 
considered here. This is a non-trivial result since relations of this kind can only be proven to hold 
in an equilibrated plasma. If generic, this implies that a kind of generalized KMS condition applies 
also in a non-equilibrium setting. 

The organization of the paper is as follows: In Section II we derive integral expressions for the 
retarded quark self-energy in a system with an anisotropic distribution obtained from contracting 
an isotropic distribution in one direction. We show plots of the different components of this self- 
energy for different anisotropy strengths and various orientations of the wave vector with respect 
to the direction of the anisotropy. We point out the strong dependence of the self-energy on the 
strength of the anisotropy and the angle of propagation with respect to the anistropy direction. 
In Section II A we prove analytically that for the case that the wave vector of the collective mode 
lies in the direction of the anisotropy there are no unstable modes. The same proof is performed 
in Section II B for the extremely anisotropic limit and arbitrary orientation of the wave vector. In 
Section III we extend our previous results to the real-time formalism and compare with the results 
obtained in the imaginary time formalism. 

II. ANISOTROPIC QUARK SELF-ENERGY 

The integral expression for the retarded quark self-energy for an anisotropic system has been 
obtained previously [40] and is given by 



Note that if there were, in fact, fermionic unstable modes one would expect extra generation of fermions 
and anti-fermions which would naively increase electromagnetic emission from the plasma. 




(1) 
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where Cp = {N^ - l)/2Nc, = /c^V(27r)^ and 

/(p) = 2(n(p)+n(p)) + 4ng(p) . 

To simplify the calculation we follow Ref. [43] and require the distribution function /(p) to be 
given by 

/(P) = fdP) = m) /iso (Vp' + e(p-fi)') • (2) 

Here n is the direction of the anisotropy, > — 1 is a parameter reflecting the strength of the 
anisotropy and N{^) is a normalization constant. For the application to heavy ion collisions fi is 
the beamline (longitudinal) direction and the relevant anisotropy parameter at times r > {pt)~^ 
is positive, ^ > 0, corresponding to an oblate distribution. 

To fix A^(^) we require that the number density to be the same both for isotropic and arbitrary 
anisotropic systems, 

/isob) = fdP) = N{0 /iso (Vp'+e(p-fi)2) , (3) 

and can be evaluated to be 

N{0 = VTK- (4) 
Using Eq. (2) and performing the change of variables 

= p^l + ■ hf) , (5) 

we obtain 

m) = mlVlTl / f (1 + ^(P • n)r' ^ , (6) 

where 

'^q = ~^J^ dPPhsoijp)- (7) 

We then decompose the self-energy into four contributions 

E(i^) = 70Eo + 7-5:. (8) 

The fcrmionic collective modes are determined by finding all four-momenta for which the de- 
terminate of the inverse propagator vanishes 

det5-^ = 0, (9) 

where 

iS-\P) = 

= l^A, . (10) 

with A{K) = {ko - Eo, k - 5]). Using the fact that det(7''A^) = (^''A^)^ and defining = a • A 
we obtain 

^0 = . (11) 
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FIG. 1: Real and imaginary part of Eq as a function of a;/fc for ^„ = 7r/4 and ^ = {0, 10, 100}. 





FIG. 2: Real and imaginary part of E^; as a function of oj/k for On = 7r/4 and ^ = {0, 10, 100}. 
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FIG. 3: Real and imaginary part of as a function of uj/k for ^„ = 7r/4 and ^ = {0, 10, 100}. 



In practice, we can define the z-axis to be in the n direction and use the azimuthal symmetry 
to restrict our consideration to the x — z plane. In this case we need only three functions instead 
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FIG. 4: Real and imaginary part of Eq as a function of u/k for ^ = 100 and On = {0, 7r/4, 7r/2}. 



of four 



R{w — k cos 6nX, ksinOnV^ — x'^) 

, Vl — x'^S(w — k cos 9nX, k sin 6nV^ — x'^) 
ax 



1 + ^x 



where 



1 ^ /•! 

^0{w,k,en,O = 2*^9 ^1+^ J ^ 

1 2 /7— T /"^ J - A; cos OnX, k sin V 1 - a;2) 

R{a, b) = a + 6 + ieV a — 6 + ze^ , 

5(a,6) = [ai?(a,&) - 1] . (13) 



In Figs. 1 through 3 wc plot the real and imaginary parts of the quark self-energies Eq, E 
and E^ for ^ = {0,10,100}. From these Figures we see that the spacelike quark self-energy is 
strongly affected by the presence of an anisotropy with a peak appearing at u/k = sin^„ for 
strong anisotropies. To further illustrate this in Fig. 4 we have plotted Eq for £^ = 100 and 

= {0, vr/4, 7r/2}. From this Figure we see that there is a large directional dependence of the 
spacelike quark self-energy. Note that this could have a measurable impact on quark-gluon plasma 
photon production during the early stages of evolution since screening of infrared divergences in 
leading order photon production amplitudes requires as input the hard-loop fermion propagator 
for spacelike momentum. Wc return to this point in Section III and sketch how to calculate 
photon emission from an anisotropic quark-gluon plasma. Assuming the necessary measurements 
of the rapidity dependence of the thermal photon spectrum could be performed, photon emission 
could provide an excellent measure of the degree of momentum-space anisotropy in the partonic 
distribution functions at early stages of a heavy-ion collision. 

For general ^ and On we have to evaluate the integrals given in Eq. (12) numerically. To find 
the collective modes we then numerically solve the fermionic dispersion relations given by Eq. (11). 
As in the isotropic case, for real timelike momenta (|a;|>|/c|, lui{uj/k) = 0) there are two stable 
quasiparticle modes which result from choosing either plus or minus in Eq. (11).^ We have looked 



Note that there are four solutions to the dispersion relations since each solution exists at both positive 
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for modes in the upper- and lower-half planes and numerically wc find none. In the next section 
we explicitly count the number of modes using complex contour integration and demonstrate that 
there are no unstable collective modes in two special cases. 



A. Special case: k || n 



Let us consider the special case where the momentum of the collective mode is in the direction 
of the anisotropy k || n, i.e., On = 0. In this case the integrals in Eq. (12) can be evaluated 
analytically. E-^ becomes zero, while the other components read 



1 

— 

2 ■"«e^2 + fc2 

S^(a;,fc,^„ = 0,0 = 1^40^^ 



2Y^a;arctan V^+fcln 



(jj + k 
uj — k 
oj + k 



-2—=k arctan \/? -|- a;ln , 



Eq. (11) simplifies to 



So = ±{k - S,) . 



(14) 
(15) 



Nyquist analysis 

We now show analytically for this special case that unstable modes do not exist. This is done 
by a Nyquist analysis of the following function: 

k,0=0J- Eo(a;, k, OT[k- ^uj, k, 0] ■ (16) 

In practice, that means that we evaluate the contour integral 

which gives the numbers of zeros minus the number of poles P of f^z in the region encircled by 
the closed path C. In Eq. (17) and in the following, we write the functions f^i in terms of 2; = uj/k 
and for clarity do not always state the explicit dependence of fzf on k and ^. Choosing the path 
depicted in Fig. 5, which excludes the logarithmic cut for real z with < 1 of the function (16), 
leads to i-* = and the left hand side of Eq. (17) equals the number of modes A'^. Evaluation of 
the respective pieces of the contour C for each /_ and /+ leads to 

iV^ = 1 + + + 1 = 2, (18) 

such that for the total number we get is = + A^_|- = 4, which corresponds to the stable modes 
(two for positive lo and two for negative lo). The four contributions in (18) are the following: 

1. The first 1 results from integration along the large circle at \z\ S> 1. 

2. The first zero is the contribution from the path connecting the large circle with the contour 
around z = ±1. 



and negative w. 
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FIG. 5: Contour C in the complex z-plane used for the Nyquist analysis. 



3. The second zero stems from the two small half-circles around z = ±1 



4. The last 1 is obtained from integration along the straight lines running infinitesimally above 
and below the cut between z = —1 and z = 1. See below for details on this integration. 



The last contribution can be evaluated using 



/l+it 
-1+ie 



dzl^ 

1+ie /=f(-2) 



In 



+ 27rin , 



(19) 



M-l + ie) 

for the line above and the corresponding expression for the line below the cut. n is the number 
of times the function /=p crosses the logarithmic cut located on the real axis, running from zero to 
minus infinity. This cut is due to the appearance of the logarithm on the right hand side of Eq. 
(19). In the sum of the line integrations above and below the cut diverging contributions from the 
first part on the right hand side of Eq. (19) cancel and we are left with a contribution of 27rz for 
each function. Furthermore it is necessary to show that neither /_ nor /+ crosses the cut. The 
proof is given in some detail for /_ and is performed analogously for /+. From Eq. (16) we find 
for /_: 



f.iz,k,0 = z-l + 



2(1+^^2) /c2 



-2 + I j A/Carctan ^/^+{z-l)\n 



z + l 
z-1 



(20) 



We want to study whether this function crosses the real axis in the range Re[z] G [—1, 1] for 
Im[2;] 0, i.e., whether the imaginary part of /_ changes sign in that range. On the straight line 
infinitesimally above the cut the imaginary part of /_ is given by 

TT VTT^ix - 1) 



Im 



e— >0 



(21) 



2 P(l + Ca;2) ' 

for real x. It is only zero for x = 1, which means that the function /_ can not cross but merely 
touch the cut within the limits of the integration. On the straight line below the cut we get the 

same result (21) with a minus sign. For /+, we find that the imaginary part in the regarded range 
only becomes zero for x = —1, which means that the logarithmic cut is not crossed within [—1,1] 
either. Hence we have proved for the case k || fi that there are no more solutions than the four 
stable modes. In particular we have shown that unstable fermionic modes can not exist. 



8 



B. Large-^ limit 



In the extremely anisotropic case where ^ ^ cx) the self-energies for arbitrary angle On can be 
calculated explicitly. The distribution function (2) becomes [63] 

lim /^(p) = S{p ■ h) r dxU ip\/l + xA . (22) 

With n in the z-direction this implies that p lies in the x-y-plane only. As in Section II, due to 
azimuthal symmetry, we consider the case where k lies in the x-z-plane only. Using (22) we obtain 
from Eqs. (12) 



2 ^ ^oj + k sin 9n\/oJ — k sin 9^ ' 

TT 2 / ^ 



2fesm6'„ ^ V + ^sm6'„Va; - A;sm6'„ / 
Since Pz is always zero, vanishes. Eq. (11) now becomes 

a; - Eo = ± - S,)2 + k^ . (24) 



Nyquist analysis 

Again, we only find four stable modes and will now show analytically that these are the only 
solutions in the large ^-limit for arbitrary angle 9n- The cut resulting from the complex square 
roots in (23) can be chosen to lie between 2; = — sin 9n and z = sin 9n on the real axis. The Nyquist 
analysis can then be performed analogously to that in Section II A with the contour in Fig. 5 
adjusted such that the inner path still runs infinitesimally close around the cut. Using this path 
in the evaluation of Eq. (17) for the functions 

k, 9n) = CO -J^o^F Vikx - S^)2 + fc2 ^ (25) 

we find the number of solutions to Eq. (24) to be 

N^ = l + + j + ^ + ^ = 2, (26) 

so that again there arc = A+ + N- = 4 solutions, which are the known stable modes. The 
decomposition in (26) is done as follows: 

1. The first contribution to N^: comes from integration along the large outer circle at \z\ S> 1. 

2. The zero stems from the paths connecting the outer and the inner circle. 

3. The two contributions of 1/4 result from integrations along the small circles around — sin^„ 
and sin 9n- 

4. The last contribution of 1/2 comes from integration along the straight lines running infinites- 
imally close above and below the cut. We discuss this part in further detail below. 

The last contribution can be obtained using Eq. (19). For the evaluation of the limit e ^ it is 
essential to note that the /qi behave like Ine or l/(lne) (depending on which function is evaluated 
on which line) and are both negative as e ^ 0. This results in a contribution of +iTr for each 
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function and integration, because in all cases the imaginary part of both functions can be shown to 
be positive in the regarded limit. All other contributions, including the diverging parts ibln(— Ine) 
cancel in the sum of the results from the upper and lower line. 

Again, we need to show that the functions /ip do not cross the logarithmic cut for z G 
[— sin0„,sin^„], i.e., that n = in Eq. (19). It is possible to find an analytic expression for 
the imaginary part of /=p using 

Imi/x + iy = sgn(y) J y/ x'^ + - x , (27) 
v2 

for the imaginary part of the square root appearing in (25) with real x and y. Then the only 
solutions to 

lm/^(z) = (28) 

are found analytically to be Re( z) — sinO^i and lie(z) — — sin^^ for f— and fj^ respectively. This 
means that the cut is not crossed during the integration along the straight lines and that the 
contribution from this piece is in fact 1/2. 



III. FERMION SELF-ENERGY FROM THE REAL-TIME FORMALISM 



In this section we extend our previous results to the real-time formalism and demonstrate that 
the high-temperature limit of the Kubo-Martin-Schwinger formula, E12 = — S21, holds even for the 
non-equilibrium configuration considered here. We will use the real-time formulation of Rcfs. [64- 
69]. In this case both propagators and self-energies become 2x2 matrices. The free propagators 
are given by 



S{K) = {t + m) 



1 






-1 

K'^ —ni^ —ie 
2 ™2 



h{K) -e{-ko) + fF{K)\ 
-e{ko) + fF{K) /f(K) / 



(29) 



with the general fermion distribution function f-p{K). 

The components (12) and (21) of the self-energy matrices are related to the emission and 
absorption probability of the particle species under consideration [67, 70, 71]. To lowest order 
photons are produced via annihilation and Compton processes 



q + q^g + -1, q{q) + g ^ q{q) + 7 ■ 
Within the real-time formalism the rate of photon emission can be expressed as [36] 

from the trace of the (12)-element II12 of the photon-polarization tensor. 



(30) 



(31) 



-e'e^AT.y ^Tr[^^iSUP)\HLl^^iS2i{P-Q) + riSi2iPh^iS*2AP-Q)\HL] , 

(32) 

where Cq is the quark charge. Here S12 and S21 are the free fermion propagators from Eq. (29) and 
propagators with an HL subscript are the full propagators in the hard-loop approximation. The 
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P(/P 

FIG. 6: 6'„-dependent part a° of the self-energy E?2. On e {0, 7r/4, 7r/2} and ^ = 100. 



hard-loop propagators satisfy a fluctuation dissipation relation, which in the quasi-static case is 
given by 

'^'12/21 (^)\hL = •S'ret(-P)l_f/L ^12/2l(-P) SI^v{P)\hL ■ (^^) 

The retarded propagator reads 

^„.(^)I..= ^_„,U(P) . (34) 

where is the retarded self-energy given in Eq. (1). The advanced propagator follows analo- 

gously with the advanced self-energy and to one loop order S12 is given by 

Si2(P) = 2ig''CF J 0^Su{K)Au{Q) , (35) 
where A12 is the (12)-element of the matrix boson propagator given by 

With the anisotropic distribution function (2) S12 can be evaluated in the hard-loop approximation 
to read 



^-(^) = ^0/^^i^^V_, ^^(1 



T 



+ Cx2)3/2 



5 (5-) NiOfFik) (iV(6/r (^) + 1) + IT ^ (9+) NiOfFik) (iV(6/r(^) - l) 

ko=k ^ y K ;i.o=-fc ^ ^ 

(37) 



where 



9± 



±Po + P (^sin On \J\ — cos -I- cos ^n^:^ j (38) 



and we chose p to lie in the x — z-plane and used the change of variables (5) for A;. Note that 
in the hard-loop limit one can ignore the quark masses and hence they have been explicitly set to 
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Q. 




FIG. 7: 6'„-dependent part a" of the self-energy Sfa- On e {0,7r/4,7r/2} and ^ = 100. For On = ^ 
= 0. 



zero above. The term k^/k does not depend on k and is given by (±l,sin6'cos^, sin6'sin^, cos^). 
Evaluation of the 5-function leads to 



with 



k^' 
T 



i 



-iS^2 (P^ 0n,O=^OC^'{P,On,O+BP'' {P, en,0, 

(1 - 



ko=k 



(1 - xfy/'^{pCOS0n+PO^Xi) - psm9nXi{l + ^) COS( 

(l-x|)V2 



feo=— fe 



(1 - xf)V2(pcos 9n - Po^Xi) - p sin 9nXi{l + ^) cos. 



where the Xj and solutions to 

k 



Po + P { sin 9nV 1 — cos ^ + cos 



-po + P [sin 9nVl — x^ cos (f) + cos 0„x^ 
= 0, respectively, and 



A 



B 



^ I dkkNiOfr{k){N{Omk) + l) , 



87r2 



J dkkNiofi'^ik) (iv(o/rw - 1) . 



(39) 

(1 - ) 

9{l-xj), 
(40) 
= and 

(41) 
(42) 



There can be A?^ € {0, 1, 2} solutions for both Xi and Xi, depending on the parameters p,po, On and 
(j). Note that ^ is also given in terms of the Xj. It is easily verified that 

a''{P,0n,C) = -f3''{P,0n,O, (43) 
such that Eq. (39) greatly simplifies to read 

-iE^2(P,e„,0 = {A-B)a''{P,0n,O, (44) 

where 



A-B = 



87r2 



/•OO 1 

iV(0 dkk{fW{k) + ff%k)) =-m^iV(e), (45) 
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FIG. 8: 6'n-dependent part of the self-energy e {0, 7r/4, 7r/2} and ^ = 100. For = 7r/2 



assuming equal quark and antiquark distributions. We did not present the analogous explicit 
calculation of S21, but find for it the same result as for ^,12 with A and B interchanged. We also 
verified that 1112 and S21 fulfill the general relation 

S2i-Si2 = 2iImS, (46) 

with the retarded self-energy S given in Sec. 11. Furthermore, since E21 is given by Eq. (44) with 
A and B interchanged it follows within the hard-loop approximation that with the form of the 
anisotropic distribution function assumed here it always holds that 

S12 = -S21 , (47) 

which can be seen as a high-tempcrature limit for the Kubo-Martin-Schwinger relation in equilib- 
rium, but also holds for finite ^ and hence for non-equilibrium. Eqs. (46) and (47) show that in 
order to calculate the hard-loop photon production rate from an anisotropic plasma one need only 
know the retarded self-energy. We plot the functions a for an anisotropy parameter of ^ = 100 
and different angles 9n in Figs. 6, 7 and 8 to emphasize the strong angular dependence once more. 



IV. CONCLUSIONS 

In this paper we have extended the exploration of the collective modes of an anisotropic quark- 
gluon plasma by studying the quark collective modes. Specifically, we derived integral expres- 
sions for the quark self-energy for arbitrary anisotropy and evaluate these numerically using the 
momentum-space rcscaling introduced in previous papers. Using direct numerical calculation we 
found only real timelike fermionic modes and no unstable modes. Additionally using complex con- 
tour integration we have proven analytically in the cases (a) when the wave vector of the collective 
mode is parallel to the anisotropy direction with arbitrary oblate anisotropy and (b) for all angles 
of propagation in the limit of an infinitely oblate anisotropy that there are no fermionic unstable 
modes. Finally, we calculated the fermion self-energy of an anisotropic plasma in the real-time 
formalism and demonstrated that within the hard-loop approximation the high-temperature limit 
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of the Kubo-Martin-Schwingcr formula, S12 = — S21, holds even for the non-equilibrium config- 
uration considered here. This means that it suffices to only have the retarded self-energy S in 
order to complete a calculation of photon production from an anisotropic plasma in the hard-loop 
framework. This calculation is currently underway [72] . 
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APPENDIX A: SMALL-,^ LIMIT 



In the limit ^ we can evaluate the quark self-energy in a power series in the anisotropy 
parameter ^. To linear order in ^ we obtain 

So = + I { I (3 cos + 1) + Sjf° [cos 20„ + 1 - (3 cos + l)z^] } , (Al) 

^,so ^ A I i (5 COS 2^„ + 3) + ZT}I° [3 cos 20^ + 3 - (5 cos 20„ + ?,)z^] | , (A2) 

I ^ (5 cos 2^„ - 1) + 3S^«° [3 cos 2^„ - 1 - (5 cos 2^„ - l)z'^] | , (A3) 



sin0,_ 

= sr + 

cos tin 

where 



cos On 12 



'm'i / ijO + k \ 
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